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Abstract 

The measurements of the statistical properties of the Cosmic Microwave Back- 
ground (CMB) fluctuations enable us to probe the physics of the very early 
Universe especially at the epoch of inflation. A particular interest lays on the 
detection of the non-Gaussianity of the CMB as it can constrain the current 
proposed models of inflation and structure formation, or possibly point out new 
models. The current approach to measure the degree of non-Gaussianity of the 
CMB is to estimate a single parameter f^if^ which is highly model-dependent. 
The bispectrum is a natural and widely studied tool for measuring the non- 
Gaussianity in a model-independent way. This paper sets the grounds for a 
full CMB bispectrum estimator based on the decomposition of the sphere onto 
projected patches. The mean bispectrum estimated this way can be calculated 
quickly and is model-independent. This approach is very flexible, allowing ex- 
clusion of some patches in the processing or consideration of just a specific region 
of the sphere. 
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1. Introduction 

Inflation is the currently favored theory of early Universe which predicts 
an early short period of rapid expansion and explains the origin of primordial 
perturbations. Many models of inflation predict weakly non-Gaussian primor- 
dial curvature perturbations and the primordial curvature perturbation ^{x) is 
parametrized using the local model [9] as follows: 

= + /^r'(*L(^) - (1) 

where f^if^ is the non-linear coupling constant in the local model. 

Due to its simplicity the local model is highly favored especially because 
all the higher-order moments are determined in terms of this f^{f^ parameter. 
Most of models of inflation only predict a value for f^ff^ to characterize the 
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non-Gaussianity of the CMB. However, there exists other models of inflation 
that predict different types of deviations from Gaussianity, detailed calculation 
of which have been investigated in [l|, [H . 

The observed CMB temperature fluctuations AT/T are related to the pri- 
mordial curvature perturbation through the following non- linear relation 

^«5t(*(x)), (2) 

where qt is the radiation transfer function. On small scales, this function is 
very complex and it is evaluated numerically by solving the Boltzmann trans- 
port equation. On very large scales, in the Sachs- Wolfe regime (/ < 10), this 
relationship simplifles to ^ = . A number of theories of inflation have been 
proposed that make different predictions about the CMB fluctuations. Measure- 
ments of the statistical properties of the CMB are a direct test of inflation which 
can help to rule out the many models of inflation that have been proposed. 

The observed CMB anisotropics AT/T can be expanded onto spherical har- 
monics: 

AT/T = Y,(^imYim- (3) 

Im 

If the CMB is Gaussian, it is fully described by its angular power spectrum: 

a = ^ E (4) 

m= — l 

But, as explained previously, many realistic models predict deviations from 
Gaussianity. Even if, the level of non-Gaussianity is predicted to be very small in 
single fleld slow-roll inflation model, there is a large class of more general models 
that predict a substantially high level of primordial non-Gaussianity. If these 
models are true, the power spectrum provides a limited insight to the physics 
of the very early Universe and higher-order estimators are needed to probe 
the CMB non-Gaussianity. The bispectrum is a natural model-independent 
approach to probe the small departure from Gaussianity that could originate 
during inflation. The CMB angular bispectrum may be calculated from product 
of three spherical harmonic coefficients of the CMB temperature fleld. For 
Gaussian flelds, the expectation value is exactly zero. Given statistical isotropy 
of the universe, the angular bispectrum Bi^i^i^ is given by (Komatsu et al, 2001): 

^hhh — 5^ ( j (Ilimi0^l2m2^l3m3^ (5) 

mi,m2,m3 ^ ^ 

where the matrix denotes the Wigner-3j symbol, Bi^i^i^ satisfy the triangle 
condition: |/i — /2I ^ ^3 ^ + ^2 and the parity invariance: h -\- h h = even. 
Then, a non-zero bispectrum is a signature of a more complicated inflationary 
period than a simple inflation slow-roll model predicts. 
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However, computation of a full bispectrum is very time-consuming and it 
is usually assumed that the signal is too weak for each of the multipoles to 
be measured individually. Instead a least squares fit to compare the observed 
bispectrum with a particular (separable) theoretical bispectrum is used. Thus, 
most non-Gaussianity studies focus on estimating the /^l^^ parameter because 
the bispectrum is fully specified by this parameter in the case of a local model 
([1]). A fast estimator for /^l^^ has been developed by [ll] and improved by [1]: 



fi 



NL 



local 

Zl</2<^3 



prim 

^ Ci^Ci^Ci^ 



where Ci is the theoretical power spectrum and Sj^™ theoretical bispec- 

trum jof for /nl^^ = 1. The statistics 5'prim defined in [8] has only a complexity 
of 0{N^^'^) whereas the full bispectrum analysis is 0{N^^'^). A detection of 
/nl^^ > 10 will rule out most of the existing inflation models. 

However, an expression like equation ([1]) is not general and there are many 
other inflationary models that predict different types of deviations from Gaus- 
sianity. Other models have been investigated and detailed calculations of other 
form of non-Gaussianity have been carried out (see for example (ll, Q). How- 
ever, all these estimators are highly model-dependent and may not be able to 
constrain, many of the non-Gaussiannity signatures. In this paper, we will in- 
vestigate a more general and model-independent approach to probe the CMB 
non-Gaussiannity by measuring the full bispectrum. If the non-Gaussianities 
are of the local type, the bispectrum will reach a maximum in case of squeezed 
configurations (i.e. one wave vector is much smaller than the other two). 

Full bispectrum estimation on spherical data using the spherical harmonics 
definition (|5]) is very time consuming. This paper tries to tackle this problem 
by introducing a promising approach for accelerating the calculation of the fully 
general bispectrum. This method is based on the projection of the sphere onto 
small-field projected maps for which a Fourier decomposition is used to esti- 
mate the bispectrum. A mean bispectrum estimator can then be obtained by 
combining results from all the patches. This paper sets the grounds for this 
method however for full optimization of the estimator, we need more accurate 
no n- Gaussian simulations which are not yet available and an analytical pre- 
diction of a full bispectrum for a given level of non-Gaussianity which is very 
complex to code up. 

The outline of the paper is as follows. In ^ the main issues of the method are 
pointed out and the method is optimized to do power spectrum estimation. In 
§2.4l the same method is used to accelerate the full bispectrum estimation. Some 
adaptations of the method are required to optimize the bispectrum estimation. 
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2. Decomposition of the sphere onto rectangular Cartesian maps 

In this section, we will describe a method to speed up the spectral analysis 
by decomposing the sphere onto patches. Such method is similar to the Welch's 
method [12, 15] that is commonly used to reduce the variance in spectral analysis 
of ID data. For a sphere, the size and the repartition of the patches have to be 
decided. Then, the pixels of the patches have to be projected onto rectangular 
Cartesian maps. A spectral estimator is then obtained by averaging the result 
of the spectral analysis on each rectangular Cartesian map. Such methods suffer 
from the limited size of the projected patches which introduces a bias. In a future 
paper, we will consider multi-taper techniques [2] that consist in averaging over 
different tapers using the full data. These methods reduce the bias since the 
data length is not shortened. 

In this section, we will follow the basic approach that consists in decomposing 
the sphere onto patches. A number of issues has to be solved before having an 
optimal estimator. 

2.1. Some issues of the method 

The first issue to be tackled is to find an optimal tessellation of the sphere. 
To facilitate the post-processing, especially the FFT required to do spectral 
analysis, the patches ought to be rectangular and preferably have the same ori- 
entation. However, no regular tessellation of the sphere by rectangles exists. 
Instead, a pseudo regular tessellation with rectangles can be obtained by allow- 
ing the patches to overlap. The overlapping also improves the power spectrum 
estimation : according to the Welch's overlapped segment method [Til, [HI, es- 
timating the power spectrum of a signal by splitting this one into overlapping 
windowed patches reduces the noise in the estimated power spectra by reduc- 
ing the frequency resolution. There still remains the size and position of the 
patches to be determined. To reduce the distortions, the size of the patches 
should depend on the map projection and is a compromise between the distor- 
tions introduced by the projection and the window function effects. 

The effect of the window function on the final power spectrum is an- 
other issue to take into account. The limited size of the patches forces us to 
use window functions. Thus, instead of analyzing the signal we will be 

analyzing the truncated signal: Sh{x^y) = s{x^y)h{x,y). In the frequency do- 
main, we obtain the following convolution product Sh{u,v) = S{u,v) * H{u^v) 
where H{u^ v) is the Fourier transform of the window function. By default the 
window function is a rectangular window, constant inside the patch and zero 
elsewhere (see the left panel of the FigE]). But, it appears that a window func- 
tion with better frequency response has to be used. The ideal window function 
is one whose frequency response is a Dirac delta function. It corresponds to an 
infinite rectangular window function which is impossible in practice. Instead, 
the frequency response normally consists of a main lobe and side lobes (see the 
right panel of the FiglS]). To be close to a Dirac function, the main lobe ought 
to be the highest and the narrowest to increase the frequency resolution and 
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the side lobes have to be the lowest to limit the mode-to-mode interaction. The 
rectangular window function is usually not recommended because of its signifi- 
cant sides lobes. 



An additional issue comes from the projection effects. The method that 
we have developed is based on the decomposition of the sphere onto a few rect- 
angular Cartesian maps. It assumes to do pixel projections from the HEALPix 
map to rectangular Cartesian maps. No matter how sophisticated the projection 
process will be, distortions are inherent in flattening the sphere. Some classes 
of map projections maintain areas, and others preserve local shapes, distances, 
and/or directions... No projection, however, can preserve all these character- 
istics. Choosing a projection thus always requires compromising accuracy in 
some way, and that is one reason why so many different map projections have 
been developed. 

Whatever the projection, if we want to keep this projection exact, we will 
have to deal with rectangular Cartesian maps having a non-regular grid. This 
will be a problem especially for the power spectrum and bispectrum estimation 
because we need to estimate the Fourier coefficients from a signal f{x^y) at 
arbitrary nodes (x^y). 

Assuming the signal / is a periodic function, it can be decomposed into 
Fourier series as follows: 



The approximation that has been used in this paper consists in finding the 
coefficients C'ni,n2(/) of ^^e Fourier series ([7]). The problem can be solved 
efficiently by noting the system matrix is a Toeplitz matrix i.e. a diagonal- 
constant matrix [see 4, for more details about the use of Toeplitz matrices in 
irregular sampling problems] 

2.2. Optimization of the method for power spectrum estimation 

As discussed in the previous section, a number of issues has to be solved in 
order to use the previous method for spectral estimation. In this section, the 
method is optimized to perform power spectrum estimation. For this purpose, 
we have generated simulations of the full sky CMB temperature as realizations 
of a random Gaussian field with a prescribed power spectrum. The cosmology 
adopted in these simulations is consistent with the WMAP parameters. We 
use the HEALPix pixelisation as WMAP and Planck missions with a resolution 
parameter of ngide = 1024. 

2.2.1. The tiling of the sphere 

Following the ideas discussed in ^ several decompositions have been tested. 




(7) 



ni = — oo n2 = — oo 
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An interesting tiling consists in decomposing the sphere onto rectangular 
patches distributed into lines of same latitude (see Fig(T]). This decomposition 




0.0 



e.o 



Figure 1: Sphere tiling using an equi- latitude decomposition 



introduces a substantial overlapping at the pole. For power spectrum estima- 
tion, this effect can be neglected by assuming the CMB field is isotropic. But 
this decomposition should not be used to detect non-Gaussianity because some 
types of non-Gaussianity can produce localized hot spots or other structures. 

Another tiling consists in decomposing the sphere onto rectangular patches 
located at the HEALPix centers of a lower resolution (see Figl2j). This decom- 
position insures a good repartition of the patches in the sky and should be 
preferred for bispectral estimation. 



Once, we have selected the optimal tiling of the sphere, the pixels of the 
patches are projected onto rectangular Cartesian maps using a projection that 
will be described in ^2.2.21 To perform the decomposition, the size of the patches 
has to be fixed in such a way that keeps the induced distortions to minimum. To 
reduce the spectral leakage, the size of the patch has to be increased to narrow 
the main lobe of the frequency response. But increasing the size of the patch 
will increase the distortions effects due to projection effects. 






Figure 2: Sphere tiling using the HEALPix centers of a lower resolution. 
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In Figini the mean power spectrum has been estimated using three different 
sizes of field: 10°xlO°, 20°x20° and 30°x30°. A Hann window has also been 
used to reduce the spectral leakage. Then, these power spectra (red) have been 
compared to the theoretical power spectrum used to simulate the full-sky CMB 
(black). 




500 1000 1 500 2000 2500 3000 



Figure 3: The effect of the size of the patches in the power spectrum estimation: 

a power spectrum (with error bars) has been estimated using three different size of patches 
(Note: a Hann window function is used to reduce the spectral leakage): 10°xlO° (top left), 
20°x20° (middle left) and 30°x30° (bottom left). Then, these three power spectra (red) have 
been compared to the theoretical power spectrum (black). The right column corresponds to 
a zoom in on the left power spectrum. 



As expected, an important spectral leakage is observed on the mean power 
spectrum estimated with patches of 10°xlO° (see top panel of Fig. [3]). With 
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patches of 30°x30° (see bottom panel of Fig. [3]), the spectral leakage is severely 
dampened but some distortions are visible at large scales, most certainly due 
to projection effects. The best result is obtained with a field of 20°x20° (see 
middle panel of Fig. [3]) which is a good compromise to both minimize the spectral 
leakage and reduce the distortions introduced by the map projection. 

2.2.2. Map projection 

The map projection introduces distortions of various classes, as it has been 
discussed in ^2.11 The choice of the map projection has to be done in such a 
way that it reduces the error in the estimated power spectrum. Two projections 
have been compared: 

1. The gnomonic projection is constructed by projecting every point of the 
sphere onto patches from the center of the sphere. Assuming the patch 
is tangent to the point S (6>o, V^o), the coordinate transformation is the 
following: 



where is the longitude and ijj is the latitude. 

2. The stereographic projection is constructed by projecting every point of 
the sphere onto patches from the sphere north pole in a plane tangent to 
the south pole. Assuming the patch is tangent to the point S (6>o, V^o), the 
coordinate transformation is the following: 



In FigHl we show the mean power spectra estimated from patches of 20°x20° 
with a Hann window using two different projections: the gnomonic projection 
(left) and the stereographic projection (right). The two mean power spectra 
(red) are compared to the theoretical power spectrum (black). The better re- 
sult is obtained with the gnomonic projection that will be used as the default 
projection for power spectrum estimation. However, we have to note that the 
map projection errors have already been reduced by fixing the size of the field 
to 20°x20° (see g22H). 

2.2.3. Windowing 

After the projection, the full-sky CMB is decomposed onto rectangular 
Cartesian maps of 20°x20°, but, as already discussed in §2.H the analysis of 
a finite signal affects the frequency analysis. The simplest way to model a patch 
of a finite size is through the usage of a rectangular window. But, this default 
window introduces an important spectral leakage. There is a lot of possible 
other window functions that can be used to reduce the spectral leakage in the 
power spectrum estimation. A few of the more common window functions have 
been compared in this study (see FiglSj): 




(8) 




(9) 
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Figure 4: The map projection effect in the power spectrum estimation: Power spectra 
estimated with two different projections: the gnomonic projection (left) and the stereographic 
projection (right). These power spectra (red) are compared to the theoretical power spectrum 
(black). The right column corresponds to a zoom of the left power spectrum. 
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Figure 5: Window functions commonly used in FFT power spectrum estimation (left) and its 
frequency response (right). A signal whose frequency is actually located at zero offset leaks 
into neighboring frequencies with the amplitude shown. The rectangular window which is 
equivalent to no windowing is the least recommended because of its large side lobes. 

1. The Rectangular window that is the default window: 

h{x,y) = l. (10) 

2. The Hann window that has been used in ^2.2.11 

h{x^y) = cos(7rx)^ cos(7ry)^. (11) 

3. The Bartlett window: 

Mx,t/) = l-|||-|||. (12) 

In Figini the mean power spectrum has been estimated by decomposing 
the CMB full-sky onto patches of 20°x20° using the three windows described 
above. These mean power spectra (in red) have been compared to the theoretical 
power spectrum that has been used to simulate the CMB full-sky (in black). 
As expected, a spectral leakage is observed with the default rectangular window 
due to the finite size of the signal (see the top panel of FigE]). The spectral 
leakage has been reduced by applying the previously described non-rectangular 
window functions (see the middle and bottom panel of Figl6]). However, the 
Hann window does the best job (see the bottom panel of Figl6]). Indeed, the 
Hann window is known to produce moderate side lobes (see the bottom panel of 
FigH]) and to have high frequency resolution which is close to an ideal window 
function (see §2.ip . Therefore, the Hann window will be used as the default 
window function for power spectrum estimation. As said previously, in this 
paper, we haven't tested the multitaper approach [2] that seems to be optimal 
to reduce the spectral leakage in the power spectrum estimated from small 
patches of the sky. This will be done in a future work. 

2.3. Validation of the method with the power spectrum 

In the previous section, we have done an optimization of the method for 
power spectrum estimation. There still room for improvements but is outside 
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Figure 6: The window function effect in the power spectrum estimation: a mean 
power spectrum (with error bars) is obtained by decomposing the sphere onto patches of 
20°x20° using a rectangular window (top left), a Bartlett window (middle left) and a Hann 
Window (bottom left). These power spectra (red) are compared to the theoretical power 
spectrum (black). The right column corresponds to a zoom in on the left power spectrum. 
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the scope of this paper. The best power spectrum estimation has been obtained 
by decomposing the sphere onto patches of 20°x20°, by projecting the pixels 
of the patch onto rectangular Cartesian maps using a gnomonic projection and 
by multiplying each projected map by a Hann window. The power spectrum 
estimated by this method is now compared with the power spectrum estimated 
from the spherical harmonics formula (equation 2]) . 




Figure 7: Mean power spectra estimated from 100 full-sky CMB simulated maps 
using two different methods: using the full-sky method based on the spherical harmonic 
coefficients (top left panel in red) and using the method based on the decomposition of the 
sphere onto patches (top right panel in pink). These two mean power spectra are compared 
to the theoretical power spectrum used to produce the simulations (in black). The bias in the 
mean power spectrum estimated with the full-sky method is about 0.5% while it is about 3% 
with the per patch method. The empirical standard deviation is given for the two methods 
in the bottom panel. 

In FiglTl we have estimated a mean power spectrum from 100 full-sky CMB 
simulated maps using the two different methods. On top left, we have the mean 
power spectrum estimated from spherical harmonics coefficients (in red) and 
in the top right, the optimized mean power spectrum obtained by the method 
based on the decomposition of the sphere onto patches (in pink). These two 
mean power spectra have been compared to the theoretical power spectrum used 
to do the simulations (in black). The two curves on the left panel lie on top of 
each other while the two curves in the right panel show a small shift. There still 
remains a small amount of leakage after windowing. This kind of bias can be 
corrected or reduced using a multitaper approach [H . The standard deviation 
is significantly smaller in the per patch method than in the full-sky method (see 
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the bottom panel of Fig. [7|), this is because the frequency resolution has been 
reduced. 

It could be thought, the per patch method is not recommended for power 
spectrum estimation because the decomposition of the sphere onto patches takes 
some time and there already exists a fast method for a fully CMB power spec- 
trum estimation. However, in practice, we never have access to a full-sky CMB 
map because of the contamination by residual foregrounds. For this reason, the 
per patch method can be preferred for some applications. 

2.4- Application of the method 

As mentioned previously, we never have access to a full-sky CMB map. The 
CMB map obtained by a method of component separation is always partially 
masked to discard contaminated pixels which introduces problems in power 
spectrum estimation. About 15-20% of the most contaminated data is removed 
mostly in the galactic plane. Several methods exist to overcome the missing 
data problem, the most common one to be the MASTER method [7] that uses 
apodization windows. The method that has been proposed in this paper is 
an alternative to methods such as MASTER that try to correct the power 
spectrum from the effect of missing data. The approach is totally different 
because the problem of missing data is solved by avoiding the patches with 
important contaminations on the power spectrum estimation. 

Following the same idea, we can imagine to use this method as a diagnosis 
of the component separation quality. An estimator of the power spectrum per 
latitude can be obtained by averaging the result of the spectral analysis on 
each patch located at the same latitude (see Fig. [T]). The contamination by 
residual foregrounds introduces distortions on the power spectrum that should 
increase close to the galactic plane. The level of contamination per latitude is 
an indicator of the efficiency of the component separation method. Different 
component separation methods can be evaluated and compared by this way 
using simulated data. 

Fig. [8] shows the effect of residual foregrounds in the power spectrum estima- 
tion. The left panel shows the power spectra estimated from a simulated CMB 
map purely Gaussian for different latitudes. The variance is only due to the 
number of patches per latitude that decreases moving to the poles (see Fig. [2]). 
This explains that the statistical variance becomes important for power spectra 
estimated for latitudes close to the poles (larger than |70°|). 

The right panel of Fig. [8] shows the power spectra per latitude estimated 
from a simulated CMB map with a significant level of residual foregrounds in 
the galactic region. This CMB map has been obtained by applying an optimized 
component separation method to simulations of Planck observations. The qual- 
ity of the component separation method can be evaluated by comparing this 
result to the previous result. As previously, we observe distortions in the power 
spectrum close to the pole due to the statistical variance (latitude larger than 
|70°|). But some important distortions are also present close to the galactic 
plane (latitude equal to 0° - red line and 9.5° - pink line) clearly due to residual 
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Figure 8: Power spectrum per latitude (top) and corresponding errors (bottom) 

(patches of 10°xlO°) estimated from a simulated Gaussian CMB map (ngide = 2048) without 
(left) and with (right) contamination of foregrounds. The north and south contributions have 
been merged. 
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foregrounds. This method can be used to select the better component separa- 
tion method and help to define the region to be masked. Obviously, this only 
could be accomplished on simulated data for which the true power spectrum is 
known. 

3. Pull bispectrum estimation 

The primary goal of this paper is to provide a method for accelerating the 
calculation of the full bispectrum that cannot be estimated with spherical har- 
monics in a reasonable time. In the previous section, we have introduced a 
promising approach that have been tested on the power spectrum estimation. 
The previous methodology will now be applied to estimate a full bispectrum. 
Obviously, the method needs now to be optimized for the bispectral estimation. 

3.1. Limitations of the study 

There will be two limitations to the optimization of this method. The first 
limitation is the resolution of the non-Gaussian simulations of CMB currently 
available. The bispectrum of a Gaussian field being null, non-Gaussian CMB 
simulations have to be used to test the validity of the method. For that purpose, 
we have used the CMB non- Gaussian simulations of local type provided at the 
following address: http://planck.mpa-garching.mpg.de/cmb/fnl-simulations, A 
non- Gaussian CMB temperature map can easily be computed with any desired 
level of non-Gaussianity (/nl^O t>y linear combination of the aim provided. More 
details about the simulations can be found in [3]. A major problem in our study 
is that the HEALPix resolution parameter of these simulations is ngide = 512 
which is quite low compared to Planck data {uside — 2048) which has limited 
our study. But, there exists no other non-Gaussian simulations of CMB publicly 
available with a better resolution. 

The second limitation is that no publicly available codes exist for a theoreti- 
cal computation of the full bispectrum for a given level of non-Gaussianity, due 
to its complexity. The comparison with theory is required to be sure that the 
processing is not introducing important errors in the mean bispectrum. An an- 
alytic prediction of the bispectrum has been given in [10] for non-Gaussianities 
of the local- type (i^^L^O but only for the equilateral configuration. 

In this paper, the comparison with the theoretical bispectrum has been 
barely done on the equilateral configuration. However, even focusing on a par- 
ticular configuration like the equilateral configuration, our method still remains 
more powerful than /^l^^ estimation methods because all modes of the equilat- 
eral bispectrum are reconstructed compared to a single parameter. 

3.2. Optimization of the method for bispectrum estimation 

In this section, the optimization of the method for bispectral estimation 
will be done by comparing the equilateral bispectrum estimated by our method 
on the CMB non-Gaussian simulations described previously with the analytic 
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prediction given in (lO|. However, the low resolution of the non- Gaussian sim- 
ulations will considerably limit this study. Consequently, more work will have 
to be done as soon as non- Gaussian simulations with a better resolution will be 
available. 

As for power spectrum estimation, a number of issues has to be solved in 
order to adapt the previous method to do bispectral estimation. The first issue 
is the tiling of the sphere. As said previously, for bispectral estimation, the 
decomposition of the sphere onto rectangular patches located at the HEALPix 
centers of a lower resolution (see Fig. [2j) is preferred because this decomposition 
ensures a good repartition of the patches in the sky. But the size of the patches 
has to be reduced significantly compared to previous decomposition because the 
bispectrum is very sensitive to the non-Gaussianities introduced by the projec- 
tion effects. As a result, the gnomonic projection will be preferred to do the 
projection from the pixels of the patches onto the rectangular Cartesian maps 
because this projection introduces very little distortions for small patches [6]. 
However, the size of the patches has to be fixed in such a way the distortions 
are reduced. In FiglHl the mean bispectrum has been estimated using three 
different size of field: 7°x7°, 10°xlO°, 17°xl7°. The pixel projection introduces 
non-Gaussianities in the projected maps that appears in the bispectrum. For 
patches of 7°x7°, the amplitude and the location of the acoustic oscillations are 
quite well detected. But, the larger the field is, the more important the projec- 
tion effects are. For patches of 17°xl7°, the CMB non-Gaussianities of the local 
type almost disappear. There is still more power at the location of the acoustic 
oscillations but the distortion effects introduce an important noise in the mean 
bispectrum. This experiment should be redone with simulations with a better 
resolution. 

Another issue is the window function. As for spectral estimation, it doesn't 
exist a perfect window function for bispectral estimation. It should be a trade-off 
between bispectral resolution and leakage effect. In Fig. [TOl we have compared 
the Hann window to the rectangular window to do bispectral estimation and 
we obtain a better result with the default rectangular window function (left). 
Indeed, the amplitude of the acoustic oscillations are better recovered with the 
default rectangular window. In a future work, a window function dedicated to 
the CMB data has to be designed to improve the bispectral estimation. Some 
authors have already tried to address the problem of finding an optimal bispec- 
tral window [see [HI, [l6| . 

3.3. Preliminary Results 

The best equilateral bispectral estimation has been obtained by decomposing 
the sphere onto patches of 7° x 7° and by projecting the pixels of the patch onto 
rectangular Cartesian maps using a gnomonic projection. The code FASTLens 
[li] publicly available has been used to compute the bispectrum in the Cartesian 
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Figure 9: The effect of the size of the patches in the bispectrum estimation: a mean 
bispectrum has been estimated using three different size of patches (7°x7°, 10°xlO°, 17°xl7°) 
from 100 full-sky non-Gaussian CMB maps (F^^l""^ = 1^0) • 




100 1 OOO 100 1000 



Figure 10: The effect of the window function in the bispectrum estimation: a mean 
bispectrum has been estimated by decomposing the sphere onto patches of 7°x7° using a 
rectangular window (left) and a Hann window (right) from 100 full-sky non-Gaussian CMB 
maps with F^^^^^ = 100 and = 500. 
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maps. This study has been hmited to the equilateral configuration of the bis- 
pectrum. As soon as a code to compute analytical prediction for the full CMB 
bispectrum will be released, the same study could be extended to all the con- 
figurations of the bispectrum. Furthermore, the study will be improved as soon 
as no n- Gaussian simulations of CMB with a better resolution will be available. 

Anyway, the equilateral bispectrum computed (see the left panel of Fig. 
[TQ|) shows the expected acoustic oscillations and the Silk damping predicted by 
theory for a CMB with non-Gaussianity of the local- type [10]. Therefore, this 
preliminary result shows that we can compute an equilateral bispectrum which 
is in good agreement with the analytical predictions. 

4. Conclusion 

The goal of this paper is to present a fast way to compute a full-sky CMB 
bispectrum. It should be pointed out that it is currently too hard to directly 
measure this full CMB bispectrum. This paper introduces a promising approach 
for accelerating the estimation of the full bispectrum on the sphere. This method 
involves the decomposition of the HEALPix map onto small projected Cartesian 
maps. A mean full bispectrum can then be estimated by combining results from 
all the projected maps. 

First, this approach has been used to estimate the power spectrum of a full- 
sky CMB map. A number of optimizations have been done to obtain the best 
power spectrum estimator. An interesting application of the method in order 
to test the quality of the CMB component separation on the galactic region has 
then been presented. 

The approach has then be applied to the full bispectrum estimation to ac- 
celerate its computation. The approach presented in this paper enables a fast 
reconstruction of the whole bispectrum directly from the observational data. 
A number of optimizations have been performed to improve the quality of the 
bispectral estimation. However, these optimizations has been only tested on the 
equilateral configuration of the bispectrum because of the lack of analytical pre- 
dictions for the full bispectrum. Anyway, this study could be easily extended to 
other configurations as soon as a code to compute analytical predictions for the 
full CMB bispectrum will be released. Another limitation of the study comes 
from the resolution of the non-Gaussian CMB simulations used for the analysis. 
This is a preliminary result that should be improved as earlier as is feasible. 

However, the equilateral bispectrum that is computed from the non- Gaussian 
CMB simulations using this approach is in very good agreement with the analyt- 
ical predictions. Indeed, the features expected by the theory are present despite 
the poor resolution of the estimated bispectrum. Thus, this approach appears 
very promising to constrain the CMB non-Gaussianity in a model-independent 
way. 
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